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The dynamics of partonic cascades may be an important 
aspect for particle production in relativistic collisions of nuclei 
at CERN SPS and BNL RHIC energies. Within the Parton- 
Cascade Model, we estimate the production of single photons 
from such cascades due to scattering of quarks and gluons 
1 9 — * 9 7i quark-antiquark annihilation q q — > j 7 
or 77, and from electromagnetic brems-strahlung of quarks 
5^57. We find that the latter QED branching process 
plays the dominant role for photon production, similarly as 
the QCD branchings q — > qg and g — > gg play a crucial role 
for parton multiplication. We conclude therefore that pho- 
tons accompanying the parton cascade evolution during the 
early stage of heavy-ion collisions shed light on the formation 
of a partonic plasma. 



Photons have remained one of the most effective probes 
of every kind of terrestrial or celestial matter over the 
ages. Thus, it is only befitting that the speculation of 
the formation of deconfined strongly interacting matter 
- some form of the notorious quark-gluon plasma - in 
relativistic heavy-ion collisions, was soon followed by a 
suggestion that it should be accompanied by charac- 
teristic production of photons. The effectiveness of pho- 
tons in probing the history of such a hot and dense matter 
stems from the fact that, after production, they leave the 
system without any further interaction and thus carry 
unscathed information about the circumstances of their 
birth. This is a very important consideration indeed, 
as the formation of a partonic plasma is likely to pro- 
ceed from a hard-scattering of initial partons, through a 
pre-equilibrium stage, to perhaps a thermally and chem- 
ically equilibrated state of hot and dense partonic mat- 
ter. In this letter we concentrate on photons coming from 
the early partonic stage in such collisions. We will not 
consider other sources of photons, e.g., those that ac- 
company the late hadronic stage, after the partons have 
hadronized. 

During the partonic stage, photons emerge from two 
different mechanisms: firstly, from collisions between par- 
tons, i.e., Compton scattering qg — ► 57 of quarks and 
gluons and annihilation qq — > 177, qq — ► 77 of quarks and 
antiquarks; secondly, from radiation of excited partons, 
i.e. electromagnetic brems-strahlung q — > 57 of time-like 
cascades initiated by quarks. Whereas the former mech- 
anism has been studied in various contexts Q , the latter 
source of photons is less explored ||, although, as we 
shall show, it is potentially much richer both in magni- 



tude and complexity. 

A dynamical description of relativistic heavy-ion col- 
lisions has been developed in the Parton Cascade Model 
(PCM) @]. The PCM is based on the parton picture 
of hadronic interactions and describes the nuclear dy- 
namics in terms of the interaction of quarks and glu- 
ons within perturbative quantum chromodynamics, em- 
bedded in the framework of relativistic quantum kinet- 
ics H . The time evolution of the system is followed by 
Monte Carlo simulation of a set of equations describing 
space-time development and renormalization-group evo- 
lution of the particles. The procedure, implemented in 
a computer code VNI [||, traces the dynamic evolution 
of partons scattering, radiating, fusing, and eventually 
clustering to pre-hadronic states that then produce the 
final-state hadron yield. VNI, the Monte Carlo imple- 
mentation of the PCM, has been adjusted on the basis 
of experimental data from e + e~ annihilation and pp (pp) 
collisions. In the past, the PCM has been used to provide 
insight B into conditions likely to be achieved at RHIC 
and LHC energies ^fs — 200 A GeV, respectively ^/s = 5 
A TeV. Only very recently, it has been found to pro- 
vide as well a reasonable description to a large body of 
particle spectra from Pb + Pb and S + S collisions at 
CERN-SPS energy ^/s ~ 18 A GeV, in contrast to the 
belief that the PCM could not be applied at energies 
v^s < 100 A GeV. A valuable advantage of the PCM is 
that it is free of assumptions about thermal or chemical 
equilibrium conditions, since the space-time evolution of 
the matter is traced causally from the moment of colli- 
sion on, so that at any point the state of the matter is 
unambiguously determined by the preceding space-time 
history. 

Prompt photons from both parton collisions and 
brems-strahlung are ideally suited to test the evolution 
of the partonic matter as described by the PCM. They 
would accompany the early hard scatterings, the ap- 
proach to thermalization and chemical equilibration if 
these are achieved at all. 

We first discuss the aspect of collisional photon produc- 
tion from parton collisions. This is straightforward, and 
is included in the PCM in terms of the elementary 2^2 
processes which yield photons, that is, the annihilation 
and Compton processes, qq —+ 57, qq —+ 77, qg — ► qj , 
the Born cross-sections of which are well-known || . In 
the PCM approach we treat these processes within per- 
turbative QCD if the momentum q± transferred in the 
parton collision is larger than some cut-off q± Q. (This 
cut-off is introduced in the cm. frame of the colliding 
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partons, and thus we may have contributions even for 
smaller transverse momentum in the nucleus-nucleus cm. 
frame.) 

Next we turn to the aspect of radiative photon emis- 
sion from parton showers initiated by collisions, which is 
a more complex issue. There are some important and in- 
teresting differences between a branching leading to pro- 
duction of photons as compared to gluons, which have 
been pointed out in detail by Sjostrand ||. 

(i) Consider an energetic quark produced by a hard scat- 
tering, which can radiate gluons and photons competi- 
tively. The branchings q — *■ qg and q — > (77 appear in the 
PCM on an equal footing and as competing processes 
with similar structures. The probability, for a quark 
to branch at some given virtuality scale Q 2 , with the 
daughter quark retaining a fraction z of the energy of 
the mother quark, is given by: 

where the first term corresponds to gluon emission and 
the second to photon emission. Thus, the relative prob- 
ability for the two processes is, 

Tq-^q-y Qem(e g ) ^ 1 , . 

VaZq'g K a s C F ~ 200 ' U 

using Oe m = 1/137, a s = 0.25, (e 2 q ) = 0.22 and C F = 
4/3. Thus, we notice that the emission of photons is 
strongly affected by the perturbative QCD effects. This 
does not mean, though, that we can simulate emission of 
photons in a QCD shower by simply replacing the strong 
coupling constant a s with the electromagnetic a em and 
the QCD color Casimir factor Cf = 4/3 by e q . One 
has to keep in mind that the gluon, thus emitted, may 
branch further either as, g — > gg or as g — > qq, implying 
that the emitted gluon has an effective non-zero mass. As 
the corresponding probability for the photon to branch 
into a quark or a lepton pair is very small, this process 
is neglected and the photon is taken to have a zero mass. 
(However, if we wish to study the dilepton production 
from the collision, this may become an important contri- 
bution g). 

(ii) The radiation of gluons from the quarks is subject 
to soft-gluon interference which is accounted for by im- 
posing an angular ordering of the emitted gluons. This is 
not so for the emitted photons. To recognize this aspect, 
consider a quark which has already radiated a number of 
hard gluons. The probability to radiate and additional 
softer gluon will get contributions from each of the ex- 
isting partons which may further branch as q — > qg or 
g —>■ gg. It is well-known |l(| that if such a soft gluon is 
radiated at a large angle with respect to all the other par- 
tons and one would add the individual contributions inco- 
herently, then the emission rate would be overestimated, 
as the interference is destructive. This happens because 
a soft gluon of a long wavelength is not able to resolve 



the individual color charges and observes only the net 
charge. The probabilistic picture of PCM is then recov- 
ered by demanding that emissions are ordered in terms of 
decreasing opening angle between the two daughter par- 
tons at each branching, i.e., restricting the phase-space 
allowed for the successive branchings. The photons, on 
the other hand, do not carry any charge and only the 
quarks radiate. Thus photons are not subject to angular 
ordering. Pictorially, the branching structure in QCD is 
'pine-tree' like, whilst in QED it is 'oak-tree' like, 
(iii) The parton emission probabilities in the QCD show- 
ers contain soft and collinear singularities, which are reg- 
ulated by introducing a cut-off scale no- This regular- 
ization procedure implies effective masses for quarks and 
gluons, 

(?) /Mo , ? (9) Mo fo s 

m °S = V T q ' m eff = Y ' ( 3 ) 

where m q is the current quark mass. Thus the gluons can- 
not branch unless their mass is more than 2m3l = fj^, 
while quarks cannot branch unless their mass is more 
than An appropriate value for /i is about 

1 GeV ; a larger value is not favored by the data, and 
a smaller value will cause the perturbative expression to 
blow up. These arguments, however, do not apply for 
photon emission, since QED perturbation theory does 
not break-down and photons are not affected by confine- 
ment forces. Thus, in principle quarks can go on emitting 
photons till their masses reduce to current quark masses 
(or, in a dense matter environment, to the correspond- 
ing 'in- medium' Debye mass). It has also been suggested 
that if the confinement forces screen the bare quarks, 
the effective cut-off can be of the order of a GeV. These 
arguments suggest that we can choose the cut-off scale 
Ho in (||) separately for the emission of photons and that 
the study of photon emission can provide valuable insight 
about confinement at work by comparing the character- 
istics of gluon versus photon radiation from quarks. 

We now turn to our results which we present in four 
examples: S + S (y/s = 20 A GeV) and Pb + Pb collisions 
(y/s ~ 18 A GeV) at CERN-SPS, as well as S + S and 
Au + Au collisions at RHIC energy yfs = 200 A GeV. 
For the following, it is useful to keep in mind, that the 
initial partonic composition of the colliding nuclei is very 
different for the two energies: at SPS, by far the main 
component are the quarks, with relative proportions N q : 
N g w 1 : 0.4, whereas at RHIC, the gluons play the 
dominant role, with N q : N g w 1 : 1.9. 

In Fig. 1 we have plotted the production of single pho- 
tons from such a partonic matter in the central rapidity 
region for Pb+Pb system at SPS energy. The dot-dashed 
histogram shows the contribution of Compton and an- 
nihilation processes mentioned above. The dashed and 
the solid histograms show the total contributions (i.e., 
including the branchings q — > (77) when the mass scale 
/Jo in (||) for the photon production is taken respectively 
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as 0.01 and 1 GeV. We see that prompt photons from 
the quark branching completely dominate the yield for 
P± < 3 GeV, whereas at larger transverse momenta the 
photons coming from the collision processes dominate. 
The reduction of uq for the q — > 57 branching is seen to 
enhance the production of photons have lower transverse 
momenta as one expects. We have also shown the pro- 
duction of single photons from pp collisions for s/s ss 24 
GeV, obtained by WA70 (H), NA24 gl), and UA6 @ 
collaborations scaled by the nuclear overlap for zero im- 
pact parameter for the collision of lead nuclei. The solid 
curve gives the perturbative QCD results for the pp 
collisions scaled similarly. The dashed curve is a direct 
extrapolation of these results to lower p±. 
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FIG. 1. The radiation of prompt photons from partonic 
matter in central collision of Pb (156 GeV/nucleon)+Pb nu- 
clei at CERN SPS. The dot-dashed histogram gives the con- 
tribution of only the collision processes. The dashed and the 
solid histograms give the contribution of the collision plus 
branchings when the fio for the q — > 57 branching is taken 
as 1 and 0.01 GeV respectively. The pp data at yfs ss 24 GeV 
scaled by the nuclear overlap function for the central collision 
of lead nuclei and the corresponding QCD predictions (solid 
curve, arbitrarily extended to lower pr) also shown for a com- 
parison, q^ denotes the transverse momentum cut-off above 
which the hard scatterings are included in the PCM. 

In Fig. 2 we have plotted the transverse momenta of 
the single photons in several rapidity bins for Pb + Pb 
and S + S systems at SPS energy. We see that the trans- 
verse spectra scale reasonably well with the ratio of the 
nuclear overlap for central collisions for the two systems 
TpbPb/Tss ~ 15.4, which is indicative of the origin of 
these photons basically from a collision mechanism. The 
slight deviation from this scaling seen at lower p± results 
in a w 20% increase in the integrated yield at central ra- 
pidities. This is a good measure of the multiple scatter- 
ings in the PCM. In fact we have found that the number 
of hard scatterings in the Pb + Pb system is w 17 times 
more than that for the S + S system, which also essen- 
tially determines the ratio of the number of the photons 
produced in the two cases. We also note that the inverse 
slope of the p± distribution decreases at larger rapidities, 



which is suggestive of the fact that the densest partonic 
system is formed at central rapidities. 

Finally in Fig. 3 we have plotted our results for S + S 
and Au + Au systems at RHIC energy in the same fash- 
ion as Fig. 2. We see that the inverse slope of the p± 
distribution is now larger and drops only marginally at 
larger rapidities, indicating that the partonic system is 
now more dense and spread over a larger range of rapid- 
ity. Even though the p± distribution of the photons is 
seen to roughly scale with the ratio of the nuclear over- 
lap functions for central collisions IauAu/Tss ~ 14.2, the 
integrated yield of photons for the Au + Au is seen to 
be only about 12 times that for the S + S system at the 
RHIC energy. We have again checked that the number 
of hard scatterings for the Au + Au system is also only 
about 12 times that for the S + S system. 
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FIG. 2. The radiation of single photons from 

S + 5(200 GeV/ A) and Pb + P6(158 GeV/A) collisions in 
different rapidity bins. The inset shows the rapidity distri- 
bution of the radiated photons. The results for the S + S 
collisions have been scaled by the ratio TpbPb/7ss ~ 15.4 
for central collisions, fio for the quark branching q — > 57 is 
taken as 0.01 GeV. 
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FIG. 3. The radiation of single photons from S + S and 
Au + Au collisions in different rapidity bins at RHIC energy 
y/s — 200 A GeV. The inset shows the rapidity distribution 
of the radiated photons. The results for the S + S collisions 
have been scaled by the ratio TauAu/Tss ~ 14.2 for central 
collisions, fio for the quark branching q — > q-y is taken as 
0.01 GeV. 

This contrasting behavior at SPS and RHIC energies 
exhibited in the figures can be understood as follows. At 
the SPS, the partonic system can reach energy densities 
up to 5 GeV/fm 3 [Q, and multiple scatterings become im- 
portant, especially for heavier colliding nuclei. At RHIC, 
the energy density can easily become twice as large, and 
the Landau-Pomeranchuk suppression may play a coun- 
teracting role. In the PCM this effect is mimicked, rather 
crudely, by inhibiting a new scattering of partons till the 
passage of their formation time after a given scattering. 
In a future publication jl5| we shall report on a possi- 
bility of seeing these competitive mechanisms at work by 
comparing results at different impact parameters for the 
same colliding nuclei, or for zero impact paramenter for 
different colliding nuclei. 

Some other observations are worthwhile commenting 
on: Firstly, recall such branchings of the partons pro- 
duced in hard collisions correspond to a next-to-leading- 
order correction in a s . These are known to be consider- 
ably enhanced for collinear emissions. The parton shower 
mechanism incorporated in the PCM amounts to includ- 
ing these enhanced contributions to all orders, instead of 
including all the terms for a given order |l6| . It may also 
be added that the first-order corrections to the Compton 
and annihilation processes in the plasma have been stud- 
ied by a number of authors |ri| ; however in the plasma, 
(q\) pa 4T 2 , thus their contribution is limited to very low 
transverse momentum. Clearly, (g^) is much larger in the 
early hard scatterings, and thus the radiations from the 
emerging partons are much more intense and also pop- 
ulate higher transverse momenta, as seen in the present 
work. The large yield of photons from the branching of 
energetic quarks preceding the formation of dense par- 
tonic matter opens an interesting possibility to look for 
a similar contribution to dilepton (virtual photon) pro- 
duction in such collisions. 

We conclude that the formation of a hot and dense 
partonic system in relativistic heavy-ion collisions may 
be preceded by a strong flash of photons following the 
early hard scatterings. Their yield will, among several 
other interesting aspects, also shed light on the extent of 
multiple scattering encountered in these collisions [ jl5| . 
However, we stress that we have only included the con- 
tribution of photons from the partonic interactions in this 
work. It is quite likely that the hadrons produced at the 
end will also interact and produce photons, as has been 
studied recently jl8| . A comparison of those results with 
the present work shows that at the SPS energy the emis- 
sion from the early hard partonic scatterings is of the 



same order as the photon production from later hadronic 
reactions, for p± < 2-3 GeV, and dominates considerably 
over the same at higher transverse momenta. 

One of us (DKS) would like to acknowledge the hos- 
pitality of Brookhaven National Laboratory, where most 
of this work was done. We also thank Dr Bikash Sinha 
for useful comments. This work was supported in part 
by the D.O.E. under contract no. DE-AC02-76H00016. 



[1] E. L. Feinberg, Nuovo Cim. 34A, 391 (1976); E. V. 
Shuryak, Phys. Lett. B78, 150 (1978); Sov. J. Nucl. Phys. 
28, 408 (1978). 

[2] See, e.g.: J. F. Owens, Rev. Mod. Phys. 59, 465 (1987); 
J. Kapusta, P. Lichard, D. Seibert, Phys. Rev. D 44, 
2774 (1991); Erratum-ibid D 47, 4171 (1993). 

[3] T. Sjostrand, Talk given at "Workshop on Photon Ra- 
diation from Quarks", Annecy, December 1991; CERN 
-TH.6369/92. 

[4] K. Geiger and B. Miiller, Nucl. Phys. B369, 600 (1992); 

K. Geiger, Phys. Rev. D 47, 133 (1993); K. Geiger, Phys. 

Rep. 258, 376 (1995); J. Ellis and K. Geiger, Phys. Rev. 

D52, 1500 (1995); Phys. Rev. D 54, 1967 (1996). 
[5] K. Geiger, Phys. Rev. D54, 949 (1996); Phys. Rev. D56, 

2665 (1997). 

[6] K. Geiger, Comp. Phys. Com. 104, 70 (1997). The latest 
versi on of the computer program VNI can be obtained 
from \ittp:/ /rhic.phys. Columbia, edu/rhic/vm , or from the 
authors. 

[7] K. Geiger and D. K. Srivastava, Phys. Rev. C56, 
2718 (1997) ; D. K. Srivastava and K. Geiger, ( 



[8] 

[9] 

[10] 

[11] 
[12] 
[13] 
[14] 

[15] 
[16] 

[17] 



[18] 



nucl- 



th/9708025| ); to appear in Phys. Lett. B; K. Geiger, ( 
QUCl- ]h/9801007| ), to appear in Nucl. Phys. A; K. Geiger 
and B. Miiller, Heavy Ion Physics 5 (1997). 
See, e.g.: E. Eichten, I. Hincliffe, K. Lane and C. Quigg, 
Rev. Mod. Phys. 56, 579 (1984). 

K. Geiger and J. I. Kapusta, Phys. Rev. Lett. 70, 1290 
(1993). 

A. Bassetto, M. Ciafaloni and G. Marchesini, Phys. Rep. 
100, 203 (1983). 

M. Bonesini et al., WA70 Coll., Z. Phys. C 38, 371 (1988). 

C. De Marzo et al., NA24 Coll. Phys. Rev. D 36, 8 (1987). 
G. Sozzi et al., UA6 Coll. Phys. Lett. B 317, 243 (1993). 
J. Cleymans, E. Quack, K. Redlich, amd D. K. Srivas- 
tava, Int. J. Mod. Phys. A10, 2941 (1995). 

D. K. Srivastava and K. Geiger, in preparation. 

QCD and Collider Physics, R. K. Ellis, W. J. Stirling, 
and B. R. Webber, Cambridge University Press (1996) 
p.157. 

R. C. Hwa and K. Kajantie, Phys. Rev. D32, 1109 (1985); 
D. K. Srivastava and B. Sinha, J. Phys. G18, 1467 (1992); 
P. K. Roy, D. Pal, S. Sarkar, D. K. Srivastava, and B. 
Sinha, Phys. Rev. C53, 2364 (1996). 
J. Cleymans, K. Redlich, and D. K. Srivastava, Phys. 
Rev. C55, 1431 (1997). 



4 



